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Abstract 

This work presents a new inclusive search for supersymmetry (SUSY) by the ATLAS experiment 
at the LHC in proton-proton collisions at a center-of-mass energy y/s = 7 TeV in final states with 
jets, missing transverse momentum and one or more isolated electrons and/or muons. The search is 
based on data from the full 201 1 data-taking period, corresponding to an integrated luminosity of 4.7 
fb _1 . Single- and multi-lepton channels are treated together in one analysis. An increase in sensitivity 
is obtained by simultaneously fitting the number of events in statistically independent signal regions, 
and the shapes of distributions within those regions. A dedicated signal region is introduced to be 
sensitive to decay cascades of SUSY particles with small mass differences ("compressed SUSY"). 
Background uncertainties are constrained by fitting to the jet multiplicity distribution in background 
control regions. Observations are consistent with Standard Model expectations, and limits are set or 
extended on a number of SUSY models. 
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I. INTRODUCTION 

Supersymmetry (SUSY) [J^ is a candidate for physics 
beyond the Standard Model (SM). If strongly interacting 
supersymmetric particles are present at the TeV scale, 
they may be copiously produced in 7 TeV proton-proton 
collisions at the Large Hadron Collider [lCj . In the min- 
imal supersymmetric extension of the Standard Model 
(MSSM) [TlT - tla ] such particles decay into jets, leptons 
and the lightest supersymmetric particle (LSP). Jets 
arise in the decays of squarks and gluinos, while leptons 
can arise in decays involving charginos or neutralinos. A 
long-lived, weakly interacting LSP will escape detection, 
leading to missing transverse momentum (p*r]J! nlss and its 
magnitude E™ lss ) in the final state. Significant E™ lss can 
also arise in scenarios where neutrinos are created some- 
where in the SUSY decay cascade. 

This paper presents a new inclusive search with the 
ATLAS detector for SUSY in final states containing jets, 
one or more isolated leptons (electrons or muons) and 
E™ lss . Previous searches in these channels have been 
conducted by both the ATLAS [11, El and CMS plUl 
collaborations. In this paper, the analysis is extended to 
4.7 fb -1 and single- and multi-lepton channels (with jets 
and E™ lss ) are treated simultaneously. A signal region 
with a soft lepton and soft jets is introduced in order 
to probe SUSY decays involving small mass differences 
between the particles in the decay chain. A new, simul- 
taneous fit to the yield in multiple signal regions and to 
the shapes of distributions within those signal regions is 
employed. Background uncertainties are constrained by 
fitting to the jet multiplicity distribution in background 
control regions. 



II. THE ATLAS DETECTOR 

The ATLAS detector [H HU consists of a tracking 
system (inner detector, ID) surrounded by a thin su- 



perconducting solenoid providing a 2 T magnetic field, 
electromagnetic and hadronic calorimeters and a muon 
spectrometer (MS). The ID consists of pixel and sili- 
con microstrip detectors, surrounded by a straw-tube 
tracker with transition radiation detection (transition ra- 
diation tracker, TRT). The electromagnetic calorimeter 
is a lead liquid-argon (LAr) detector. Hadronic calorimc- 
try is based on two different detector technologies, with 
scintillator-tiles or LAr as active media, and with ei- 
ther steel, copper, or tungsten as the absorber mate- 
rial. The MS is based on three large superconducting 
toroid systems arranged with an eight-fold azimuthal coil 
symmetry around the calorimeters, and three stations of 
chambers for the trigger and for precise position mea- 
surements. The nominal pp interaction point at the cen- 
ter of the detector is defined as the origin of a right- 
handed coordinate system. The positive x-axis is defined 
by the direction from the interaction point to the center 
of the LHC ring, with the positive y-axis pointing up- 
wards, while the beam direction defines the z-axis. The 
azimuthal angle <f> is measured around the beam axis and 
the polar angle 9 is the angle from the z-axis. The pseu- 
dorapidity is defined as r\ = — lntan(0/2). Transverse 
coordinates, such as the transverse momentum, p^, are 
defined in the (x-y) plane. 

III. SUSY SIGNAL MODELING AND 
SIMULATED EVENT SAMPLES 

The SUSY models considered are MSUGRA/CMSSM 
[H IH, minimal GMSB |3^fioj and a number of sim- 
plified models (Mill- The MSUGRA/CMSSM model 
is characterized by five parameters: the universal scalar 
and gaugino mass parameters mo and m.1/2, a universal 
trilinear coupling parameter Ao , the ratio of the vacuum 
expectation values of the two Higgs doublets tan/3, and 
the sign of the Higgsino mass parameter [i. In this analy- 
sis, the values of mo and m 1 / 2 are scanned, and the other 
parameters are fixed as follows: tan/3 = 10, Aq = and 
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TABLE I. Simulated background event samples used in this analysis, with the corresponding production cross sections. The 
notation LOxK indicates that the process is calculated at leading-order and corrected by a factor derived from the ratio of 
NLO to LO cross sections for a closely related process. The tt, W+ light-jets and Z+ light-jets samples are normalized using 
the inclusive cross sections; the values shown for the W+ light-jets and Z+ light-jets samples are for a single lepton flavor. 
The single-top cross sections are listed for a single lepton flavor in the s- and i-channels. Further details are given in the text. 



\x is taken to be positive. A diagram showing the decay 
of the associated production of squark and gluino is de- 
picted in Fig. [T] (a). Other diagrams representative for 
the SUSY models discussed in the following are shown in 
Fig.rTJ(b-d). 

The minimal GMSB model has six parameters: the 
SUSY breaking scale A, the mass scale of the messen- 
ger fields M mcs , the number of messenger fields N 5 , the 
scale of the gravitino coupling C grav , the ratio of the vac- 
uum expectation values of the two Higgs doublets tan/3, 
and the sign of the Higgsino mass parameter /i. For the 
minimal GMSB model, the parameters tan /3 and A are 
scanned and the other parameters are assigned fixed val- 
ues: M mcs = 250 TeV, N 5 = 3, C grav = 1 and the sign 
of /i is taken to be positive. The mass scale of the col- 
ored superpartners is set by the parameter A, while the 
next-to-lightest SUSY particle (NLSP) is determined by 
a combination of A and tan (3. At low values of A, the 
NLSP is the lightest neutralino (xi) while at the higher 
values of A where this search provides new sensitivity, 
the NLSP is a stau for tan/3 > 10 and a slepton of the 
first and second generation otherwise. The NLSP de- 
cays into its SM partner and a nearly massless gravitino. 
The gaugino and sfermion masses are proportional to N$ 
and y/Ns, respectively. The parameter C grav determines 
the NLSP lifetime, set here such that all NLSPs decay 
promptly. 

Several simplified models are considered in this pa- 
per. In the "one-step" models, SUSY production pro- 



ceeds via either pp — > gg or pp — > qhq*^ where only 
left-handed squarks of the first- and second-generation 
are considered. The gluino decays to the neutralino 
LSP via g —> qq'xt ~^ qq'W ± Xii an d the squark via 
ql —> q'xt q'W^Xi, where the W-boson can be real 
or virtual. The gluino and LSP masses are varied while 
the chargino mass is set to be halfway between them. In 
a variant of the one-step model, the LSP mass is held 
fixed at 60 GeV while the gluino (squark) and chargino 
masses are scanned. 

In the "two-step" models, SUSY production proceeds 
via either pp — > gg or pp — > <7l<?2' again where squarks 
of the first- and second-generation are considered. In the 
first class of two-step models all squarks and gluinos de- 
cay via a chargino: g — > qq'xt an d Ql q'xt- The 
charginos decay via X\ ~^ or Xi ~> u ^Li m case 
of third generation sleptons, the decay to the stau is 
via Xi - > v t\- All three generations of sleptons and 
sneutrinos are allowed with equal probability, resulting 
in equal branching ratio to sleptons and to sneutrinos. In 
the second class of two-step models, the gluinos or left- 
handed squarks decay either via a chargino (g — > qq'xt or 
qh —> q'xt) 01 v i a a neutralino (g — >• qqX2 or QL — > qx_2 )• 
The events are generated such that one chargino and one 
neutralino are always present in the decays of the pair 
produced gluinos or left-handed squarks. Neutralino de- 
cays proceed via either x° &l or X2 ^ As in the 
first two-step model, all three generations of sleptons and 
sneutrinos are allowed with equal probability, resulting 
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FIG. 1. Representative diagrams for the different SUSY models considered in this analysis: (a) MSUGRA/CMSSM model 
with pp — > qg and subsequent decay of the squark via a chargino; (b) GMSB model with pp — > qq* and subsequent decay via 
sleptons and staus; (c) one-step simplified model with pp — > qLql an d subsequent decay via charginos; (d) two-step simplified 
model with pp — >■ gg and subsequent decays via charginos and sleptons or sneutrinos. 



in a 50% branching ratio to sleptons and to sneutrinos. 
Finally, in the third class of two-step models without 
intermediate sleptons, the gluino and squark decay via 
g — » qq'xf or Ql — ► q'xf'i the decay of the chargino 



then proceeds via xf -> W W± X2 ~> W^Z^Xi- This 
signature is realized in the MSSM in a parameter region 
where additional decay modes, not contained in the sim- 
plified model, may lead to a significant reduction of the 
cross section times branching fraction of the W Z signa- 
ture. 

In the first two types of two-step models, the chargino 
and neutralino have equal masses (again set to be halfway 
between the gluino/squark and LSP mass); the slepton 
and sneutrino masses are set to be equal and halfway 
between the chargino/neutralino and LSP masses. In 
the third two-step model, the xt mass is set halfway 
between the gluino/squark and LSP while the x.2 m ass 
is set halfway between the chargino and LSP. In all the 
simplified models, the superpartners that have not been 
mentioned are decoupled by setting their masses to multi- 
TeV values. 

Simulated event samples are used for estimating the 
signal acceptance, the detector efficiency, and for estimat- 
ing many of the backgrounds (in most cases in association 
with data-driven techniques). The MSUGRA/CMSSM 
and minimal GMSB signal samples are generated with 
Herwig+-|- 2.5.2 @ and MRST2007LO* Qparton dis- 
tribution functions (PDFs); ISAJET 7.80 g| is used 
to generate the physical particle masses. The simplified 
models are generated with one extra jet in the matrix el- 
ement using MADGRAPH5jp, interfaced to PYTHIA 
[11, with the CTEQ6L1 [47j] PDF set; MLM match- 
ing 48| is done with a scale parameter that is set to 
one-fourth of the mass of the lightest sparticle in the 
hard-scattering matrix element. Signal cross sections are 
calculated in the MSSM at next-to-leading order in the 
strong coupling constant, including the resummation of 
soft gluon emission at next-to-leading-logarithmic accu- 
racy (NLO+NLL) dUm. 

The simulated event samples for the SM backgrounds 
are summarized in Table HI The ALP GEN and MAD- 
GRAPH samples are produced with the MLM matching 
scheme. The ALPGEN samples are generated with a 



number of partons < N pax t on < 5 in the matrix ele- 
ment, except for W + light-flavored jets which are gen- 
erated with up to 6 partons. The Wbb : Wcc and Wc 
cross sections shown are the leading-order values from 
ALPGEN multiplied by a K-factor of 1.2, based on the 
K-factor for light-flavored jets. For the final result, mea- 
sured cross sections are used for the W/Z+ heavy-flavor- 
jets samples • The overlap between the heavy- flavored 
and light-flavored W/Z+jets samples is removed. The 
cross section for Z+jets with 10 GeV < m« < 40 GeV is 
obtained by assuming the same K-factor as for mu > 
40 GeV. The single-top cross sections are taken from 
MC@NLO; for the s- and i-channels, they are listed for 
a single lepton flavor. 

The theoretical cross sections for W+jets and 
Z+jets are calculated with FEWZ [H] with the 
MSTW2008NNLOj5| PDF set. For the diboson cross 
sections, MCFM [jjflfwith the MSTW2008NLO PDFs is 
used. The ti cross section is calculated with HATHOR 
1.2 [H| using MSTW2008NNLO PDFs. The tt+W cross 
section is taken from Ref. [32|. The ti+Z cross section is 
the leading-order value multiplied by a K-factor deduced 
from the NLO calculation at -/I = 14 TeV [33T ]. 

Parton shower and fragmentation processes are sim- 
ulated for the ALPGEN and MC@NLO samples using 
HERWIG |i| with JIMMY [H^ for underlying event 
modeling; PYTHIA is used for the AcerMC single-top 
sample and ti+W/Z. The PDFs used in this analysis are: 
CTEQ6L1 for the ALPGEN and MADGRAPH samples, 
CT10 [E3 for MC@NLO, and MRSTMCal (LO**) [H] 
for HERWIG. The underlying event tunes are the AT- 
LAS AUET2B.LO** tunes mt 

The detector simulation 60] is performed using 
GEANT4 [6l|. All samples are produced with a range 
of simulated minimum-bias interactions overlaid on the 
hard-scattering event to account for multiple pp interac- 
tions in the same beam crossing (pile- up). The overlay 
also treats the impact of pile-up from beam crossings 
other than the one in which the event occurred. Correc- 
tions are applied to the simulated samples to account for 
differences between data and simulation for the lepton 
trigger and reconstruction efficiencies, momentum scale 
and resolution, and for the efficiency and mis-tag rates 
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for 6-quark tagging. 



jets which overlap with pre-selected electrons within a 
distance AR = (Aij) 2 + (Acj)) 2 = 0.2 are discarded. 



IV. OBJECT RECONSTRUCTION 

This analysis is based on three broad classes of event 
selection: i) a hard single- lepton channel that is an exten- 
sion to higher masses of the previous search [16[ , ii) a soft 
single-lepton channel geared towards SUSY models with 
small mass differences in the decay cascade, and Hi) a 
multi-lepton channel aimed at decay chains with higher 
lepton multiplicities. The event selection requirements 
are described in detail in Sec. I VII Here the final-state 
object reconstruction and selection are discussed. 



A. Object Preselection 

The primary vertex [62j | is required to be consistent 
with the beam spot envelope and to have at least five 
associated tracks; when more than one such vertex is 
found, the vertex with the largest summed |pt| 2 of the 
associated tracks is chosen. 

Electrons are reconstructed from energy clusters in 
the electromagnetic calorimeter matched to a track in 
the ID [63[ . Pre-selected electrons are required to have 
\r)\ < 2.47 and pass a variant of the "medium" selec- 
tion defined in Ref. that differs mainly in having 
a tighter track-cluster matching in 77, stricter pixel hit 
requirements, additional requirements in the TRT, and 
tighter shower-shape requirements for |rj| > 2.0. These 
requirements provide background rejection close to the 
"tight" selection of Ref. [63[ with only a few percent 
loss in efficiency with respect to "medium" . Pre-selected 
electrons are further required to pass a pt requirement 
depending on the analysis channel: 10 GeV for the hard- 
lepton and multi-lepton channels, and 7 GeV in the soft- 
lepton channel. 

Muons are identified either as a combined track in the 
MS and ID systems, or as an ID track matched with a 
MS segment [6J, |65[ . Requirements on the quality of the 
ID track are identical to those in Ref. Pre-selected 
muons are required to have |?7| < 2.4 and a pt require- 
ment that depends on the analysis channel: 10 GeV for 
the hard-lepton and multi-lepton channels, and 6 GeV in 
the soft-lepton channel. 

Jets are reconstructed using the anti-fct algorithm 
[o6l with a radius parameter R — 0.4. Jets arising 
from detector noise, cosmic rays or other non-collision 
sources are rejected [68|. To account for the differ- 
ences between the calorimeter response to electrons and 
hadrons, pt- and ^-dependent factors, derived from sim- 
ulated events and validated with test beam and collision 
data, are applied to each jet to provide an average energy 
scale correction [68| back to particle level. Pre-selected 
jets are required to have pr > 20 GeV and \r]\ < 4.5. 
Since electrons are also reconstructed as jets, pre-selected 



B. Signal Object Selection 

For the final selection of signal events, "signal" elec- 
trons are required to pass a variant of the "tight" se- 
lection of Ref. [631 ] . providing 1-2% gain in efficiency 
and slightly better background rejection. Signal elec- 
trons must have | ^7 1 < 2.47 and a distance to the closest 
jet AR > 0.4. They are also required to satisfy isola- 
tion criteria: the scalar sum of the pt of tracks within a 
cone of radius AR = 0.2 around the electron (excluding 
the electron itself) is required to be less than 10% of the 
electron p^ ■ 

Muons in the final selection ("signal" muons) are re- 
quired to have |?7| < 2.4 and AR > 0.4 with respect to 
the closest jet. Further isolation criteria are imposed: 
the scalar sum of the pt of tracks within a cone of ra- 
dius AR = 0.2 around the muon candidate (excluding 
the muon itself) is required to be less than 1.8 GeV. The 
Pt requirements for signal electrons and muons depend 
on the signal regions and are described in Sec. IVI1 

Signal jets are required to have pr > 25 GeV and |?7| < 
2.5. In addition, they are required to be associated with 
the hard-scattering process, by demanding that at least 
75% of the scalar sum of the pt of all tracks associated 
with the jet come from tracks associated with the primary 
vertex of the event. Jets with no associated tracks are 
rejected. The above requirements are applied to cope 
with the high pile- up conditions of the 2011 data-taking, 
in particular the later part of the run. 

The missing transverse momentum is computed as the 
negative of the vector sum of the pt of all pre-selected 
electrons, pre-selected muons and pre-selected jets (after 
removing those overlapping with pre-selected electrons), 
and all calorimeter clusters with < 4.9 that are not 
associated with any of the above-mentioned objects. 

For approximately 20% of the 2011 data-taking pe- 
riod, an electronics failure created a region in the elec- 
tromagnetic calorimeter, located at < rj < 1.4 and 
—0.8 < <j) < —0.6, where no signals could be read out. 
Events with an electron in this region are vetoed for the 
entire dataset, leading to an acceptance loss of less than 
1% for signal events in the signal region. For jets, the 
amount of transverse energy (£/r) lost in the dead region 
can be estimated from the energy deposited in the neigh- 
boring calorimeter cells. If this lost Et projected along 
the £"™ lss direction amounts to more than 10 GeV and 
constitutes more than 10% of the E™ lss , the event is re- 
jected. The effect of the electronics failure is described in 
the detector simulation, and the loss of signal acceptance 
from this requirement is negligible. 

Jets arising from 6-quarks are identified using informa- 
tion about track impact parameters and reconstructed 
secondary vertices [691 ]: the 6-tagging algorithm is based 
on a neural network using the output weights of the Jet- 
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Fitter+IP3D, IP3D, and SV1 algorithms (denned in Ref. 
[69l p as input. The 6-tagging requirements are set at an 
operating point corresponding to an average efficiency 
of 60% for 6-jets in simulated tt events, for which the 
algorithm provides a rejection factor of approximately 
200-400 for light-quark and gluon jets (depending on the 
px of the jet) and a rejection of approximately 7-10 for 
charm jets. 



V. TRIGGER AND DATA COLLECTION 

The data used in this analysis were collected from 
March through October 2011, during which the instanta- 
neous luminosity of the LHC reached 3.65 x 10 33 cm _2 s -1 . 
The average number of interactions per beam crossing 
ranged from approximately 4 to 16 during the run, with 
an average of 10. After the application of beam, detec- 
tor, and data-quality requirements, the total integrated 
luminosity is 4.7 fb _1 . The uncertainty on the luminosity 
is determined to be 3.9% [zSlnl. 

Three types of triggers were used to collect the data: 
electron, muon and E^ lss . The electron trigger se- 
lects events containing one or more electron candidates, 
based on the presence of a cluster in the electromagnetic 
calorimeter, with a shower shape consistent with that of 
an electron. The transverse energy threshold at the trig- 
ger level was either 20 GeV or 22 GeV, depending on the 
instantaneous luminosity. For signal electrons satisfying 
px > 25 GeV, the trigger efficiency is in the plateau region 
and ranges between 95% and 97%. In order to recover 
some of the efficiency for high-px electrons during run- 
ning periods with the highest instantaneous luminosities, 
events were also collected with an electron trigger with 
looser shower shape requirements but with a px threshold 
of 45 GeV. 

The muon trigger selects events containing one or more 
muon candidates based on tracks identified in the MS and 
ID. The muon trigger px threshold was 18 GeV. During 
running periods with the highest instantaneous luminosi- 
ties, the trigger requirements on the number of MS hits 
were tightened; in order to recover some of the resulting 
loss in efficiency, events were also collected with a muon 
trigger that maintained the looser requirement on the 
number of hits chambers but that required in addition 
a jet with px greater than 10 GeV. This jet requirement 
is fully efficient for jets with offline calibrated px greater 
than approximately 50 GeV. The muon triggers reach 
their efficiency plateaus below a signal muon px thresh- 
old of 20 GeV. The plateau efficiency ranges from about 
70% for \r]\ < 1.05 to 88% for 1.05 < \r)\ < 2.4. 

The E™ as trigger bases the bulk of its rejection on 
the vector sum of transverse energies deposited in pro- 
jective trigger towers (each with a size of approximately 
A77 x A(f> ~ 0.1 x 0.1 for I77I < 2.5 and larger and less 
regular in the more forward regions). A more refined cal- 
culation based on the vector sum of all calorimeter cells 
above threshold is made at a later stage in the trigger 



processing. The trigger required E™ lss > 60 GeV, reach- 
ing its efficiency plateau for offline calibrated E™ lss > 
180 GeV. The efficiency on the plateau is close to 100%. 

VI. EVENT SELECTION 

Two variables, derived from the kinematic properties 
of the reconstructed objects, are used in the event se- 
lection. The transverse mass (tox) computed from the 
momentum of the lepton (£) and the missing transverse 
momentum (p^ 1183 ), defined as 

to t = \/2p e T E™ iss (l - cos(A0(f,p T miss ))), 

is useful in rejecting events containing a single W boson. 
The inclusive effective mass (m^g) is the scalar sum of 
the pt of the leptons, the jets and E™ lss : 

N isp N jet 

i=i j=i 

where the index i runs over all the signal leptons and j 
runs over all the signal jets in the event. The inclusive 
effective mass is correlated with the overall mass scale of 
the hard-scattering process and provides good discrim- 
ination against the SM background, without being too 
sensitive to the details of the SUSY decay cascade. The 
analysis in Ref. [16[ used the three or four leading-px jets 
in the calculation of the effective mass; the additional jets 
used here improve the discrimination between signal and 
background. A second definition for the effective mass, 
denoted by m c g, is based on the sum over the 2-, 3-, or 
4-leading px jets, depending on the minimum number of 
jets required in a given signal region. This variable is 
used to compute the ratio £ , ™ lss /m e ff which reflects the 
fluctuations in the E™ lss as a function of the calorimeter 
activity in the event; the definition used here improves 
the rejection of the background from mismeasured jets. 

This analysis is based on five signal regions, each tai- 
lored to maximize the sensitivity to different SUSY event 
topologies: 1,2) Signal regions requiring a hard lepton 
plus 3- or 4-jets are extensions of the previous analy- 
sis [l6( to higher SUSY mass scales; these signal regions 
have been optimized for the MSUGRA/CMSSM model 
as well as for the bulk of the one-step simplified models 
with large mass difference (Am) between the gluino and 
the LSP; 3) A soft-lepton signal region targets the sim- 
plified models with small Am, where the hard leading 
jet comes from initial-state radiation (ISR); 4) A multi- 
lepton signal region with > 2 jets is tailored to GMSB 
models; 5) A multi-lepton signal region with > 4 jets is 
geared towards the two-step simplified models with in- 
termediate sleptons and sneutrinos. These signal regions 
are described in more detail and summarized in Table [TT1 

1. Hard lepton plus three jets. Events are selected 
with the electron and muon triggers. The num- 
ber of signal leptons with px > 25 (20) GeV for 
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> 10 


> 10 




> 3 


> 4 


> 2 


> 2 


> 4 


P J T ct (GeV) 


> 100, 25, 25 


> 80, 80, 80, 80 


> 130,25 


900 900 


■> 50 50 50 50 


p add.jet (QeV) 


< 80 


— 


— 


< 50 




E^ iss (GeV) 


> 250 


> 250 


> 250 


> 300 


> 100 


m T (GeV) 


> 100 


> 100 


> 100 








> 0.3 


> 0.2 


> 0.3 




0.2 


m% c (GeV) 


> 1200 


> 800 






> 650 



TABLE II. Overview of the selection criteria for the signal regions used in this analysis. The pt selections for leptons are given 
for electrons (muons). 



electrons (muons) is required to be exactly one. 
Events containing additional signal leptons with 
Pt > 10 GeV are rejected. The number of sig- 
nal jets is required to be > 3, with a leading jet 
satisfying pt > 100 GeV and the other jets hav- 
ing pt > 25 GeV. Events with four or more jets 
are rejected if the fourth jet has pt > 80 GeV; this 
requirement keeps this signal region disjoint from 
the 4-jet signal region. In addition, the following 
conditions are imposed: tot > 100 GeV, E^ 133 > 
250 GeV, E™ i33 /m cS > 0.3, and to" c > 1200 GeV. 

2. Hard lepton plus four jets. The lepton requirements 
are the same as in the previous signal region. The 
number of signal jets is required to be > 4, with 
the four leading jets satisfying pt > 80 GeV. In 
addition, the following requirements are applied: 
to t > 100 GeV, E™ 33 > 250 GeV, E™ iss /m cS > 
0.2, and to^ c > 800 GeV. 

3. Soft-lepton selection. Events are selected with the 
i?™ lss trigger. The number of signal leptons (elec- 
tron or muon) is required to be exactly one. Elec- 
trons are required to have 7 GeV < pt < 25 GeV, 
and muons are required to be in the range 6 GeV < 
Pt < 20 GeV. Events containing an additional sig- 
nal electron (muon) with pt > 7 (6) GeV are re- 
jected. The number of signal jets is required to 
be > 2, with the leading jet satisfying pt > 130 
GeV and the second jet having pt > 25 GeV. In ad- 
dition, the following conditions are required: tot > 
100 GeV, E™ 33 > 250 GeV, and E^/m^ > 0.3. 
No explicit requirement on m"g is applied. 

4. Multi-lepton plus two jets. Events are selected with 
the electron and muon triggers. Two or more signal 
leptons are required, with a leading electron (muon) 
with pt > 25 (20) GeV and sub-leading leptons 



with pt > 10 GeV. The two leading leptons must 
have opposite charge. At least two signal jets with 
Pt > 200 GeV are required. Events with four or 
more signal jets are rejected if the fourth leading jet 
has pt > 50 GeV; this requirement keeps this sig- 
nal region disjoint from the multi-lepton plus 4-jet 
signal region. In addition the E™ 133 is required to 
be > 300 GeV. No explicit requirements are made 
on E™ lss /m e & or m 1 ^ . 

5. Multi-lepton plus four jets. The lepton require- 
ments are the same as in the multi-lepton plus two 
jets signal region. At least four signal jets with 
Pt > 50 GeV are required. In addition, the follow- 
ing requirements are imposed: E T llss > 100 GeV, 
E^ iss /m cS > 0.2, and rr^f > 650 GeV. 

In contrast to the previous analysis [l6j , no require- 
ment on the azimuthal angle between the E^ 3 vector 
and any of the jets is imposed as the background from 
multijet events is already low. This adds sensitivity to 
SUSY decay chains where the LSP is boosted along the 
jet direction. 



VII. BACKGROUND ESTIMATION 

The dominant sources of background in the single- 
lepton channels are the production of semi- and fully- 
leptonic tt events, and lU+jets where the W decays lep- 
tonically. For the multi-lepton channels, the main back- 
ground sources are Z+jets and tt. Other background 
processes which are considered are multijets, single-top, 
dibosons and tt plus vector boson. 

The major backgrounds are estimated by isolating each 
of them in a dedicated control region, normalizing the 
simulation to data in that control region, and then using 
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the simulation to extrapolate the background expecta- 
tions into the signal region. The multijet background is 
determined from the data by a matrix method described 
below. All other (smaller) backgrounds are estimated en- 
tirely from the simulation, using the most accurate the- 
oretical cross sections available (Table H}. To account for 
the cross-contamination of physics processes across con- 
trol regions, the final estimate of the background is ob- 
tained with a simultaneous, combined fit to all control 
regions, as described in Sec. HXl 

Several correction factors are applied to the simula- 
tion. The pt of the Z boson is reweighted based on a 
comparison of data with simulation in an event sample 
enriched in Z+jets events. The same correction factor is 
applied to W boson production and improves the agree- 
ment between data and simulation in the E™ lss distri- 
bution. Other correction factors are derived during the 
combined fit. The relative normalization of the ALPGEN 
samples (IV+jets, Z+jets and tt) with different numbers 
of partons (-/Vp art0 n) in the matrix element is adjusted by 
comparing the jet multiplicity distributions in data and 
simulation in all control regions. A common set of cor- 
rections is obtained for the IV+jets and Z+jets samples, 
and a separate set of common corrections is obtained for 
semi-leptonic and fully-leptonic tt decays. Neither the 
reweighting based on the px distribution of the Z boson 
nor the Ar par ton weights are applied in Figs. [2E1 below. 



A. TV/Z+jets and tt Control Regions 

The IV+jets and tt processes are isolated in control 
regions defined by the following requirements. For the 
hard single-lepton channel, > 3 jets are required, with 
a leading jet px > 80 GeV and the other jets above 25 
GeV. The lepton requirements are the same as in the sig- 
nal region. The E^ lss is required to be between 40 and 
150 GeV while the transverse mass is required to be be- 
tween 40 and 80 GeV. Furthermore, the m 1 ^ requirement 
is relaxed to be > 500 GeV. The IV+jets and tt control 
regions are distinguished by requirements on the number 
of 6-tagged jets. For the TV+jets control region, events 
are rejected if any of the three highest px jets is 6-tagged; 
the rejected events then define the tt control region. Ta- 
ble IIIII summarizes the control region definitions; Fig. [2] 
shows the composition of the TV+jets and tt control re- 
gions as a function of m^g and of the jet multiplicity A 
discrepancy between simulation and data can be seen in 
the mj,g c distribution and is discussed in Sec. IVII Bl 

For the soft-lepton channel, the control region require- 
ments on the leptons and jets are the same as in the 
signal region. However, the E!f lss is required to be be- 
tween 180 GeV and 250 GeV and the transverse mass to 
be between 40 GeV and 80 GeV. The tighter £f liss re- 
quirement, compared to the hard single-lepton control 
regions, is dictated by the trigger selection for this chan- 
nel. Again, the IV+jets and tt control regions are dis- 
tinguished by the presence of ^-tagged jets. For W+jets, 



events are rejected if any of the two highest px jets is 
6-tagged; the rejected events form the tt control region. 
Figure[3]shows the composition of the IV+jets and tt con- 
trol regions for the soft-lepton channel as a function of 
E™ lsa /m c g and the jet multiplicity. 

For the multi-lepton channels, the Z+jets control re- 
gion is defined by requiring > 2 jets with the two lead- 
ing jets having px > 80 GeV and 50 GeV, respectively, 
or with four leading jets having px > 50 GeV. In addi- 
tion, E™ 1SS < 50 GeV and an opposite-sign, same-flavor 
dilcpton pair with invariant mass between 81 GeV and 
101 GeV are required. The lepton selection requirements 
are the same as in the signal region. The tt control re- 
gion is defined with the same jet requirements as the 
Z+jets control region; at least one jet is required to 
be 6-tagged. In addition, E™ lss between 30 GeV and 
80 GeV and a dilepton invariant mass outside the win- 
dow [81,101] GeV are required. Figure 2] (top) shows the 
composition of the Z+jets and tt control regions for the 
multi-lepton channel as a function of m 1 ^ . 



B. Reweighting of W+jets and Z+jets Simulated 
Samples 

The samples of simulated IV+jets and Z+jets events 
are reweighted as a function of the generated px of the 
vector boson. A common set of corrections to the px of 
the vector boson, applied to both W+jets and Z+jets 
samples, is found to improve the agreement between data 
and simulation for a number of variables (E™ lss , m 1 ^, 
and jet px). 

The p^ distribution is measured in data by selecting a 
sample with two oppositely-charged, same-flavor leptons 
with an invariant mass between 80 GeV and 100 GeV, 
> 3 signal jets with p T >25 GeV, and m^f > 400 GeV. 
The p^ distribution in five bins of reconstructed px is 
compared to the ALPGEN simulation in five bins of gen- 
erated px, with the first four bins ranging from to 200 
GeV and the last bin integrated above 200 GeV; the ratio 
of the two distributions is taken as the p^' gGn -dependent 
weighting factor. The simulation employed here uses the 
cross sections listed in Table |TJ Only the systematic un- 
certainty from the jet energy scale is considered (in ad- 
dition to statistical uncertainties) when computing the 
uncertainty on the weighting factors. 

Figure (top) shows the pfi distribution before the 
application of the reweighting factors and after the final 
fit to all background control regions (described in Sec. 
lIXp . which includes the reweighting. The bottom half of 
the figure shows the E™ lss distribution in the hard-lepton 
W+jets control region (with the lower requirement on 
E™ lss set to 50 GeV and the upper requirement removed) . 
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multi-lepton 




W CR 
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> 2 
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> 130,25 
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> ^0 ^0 ^0 ^0 


> ^0 ^0 ^0 ^0 


TVjct (6-tagged) 





> 1 





> 1 




> i 


E™ iss (GeV) 


[40,150] 


[40,150] 


[180,250] 


[180,250] 


< 50 


[30,80] 


m T (GeV) 


[40,80] 


[40,80] 


[40,80] 


[40,80] 






m% c (GeV) 


> 500 


> 500 










m u (GeV) 










[81,101] 


< 81 or > 101 



TABLE III. Overview of the selection criteria for the W+jets, Z+jets and tt control regions (CR). Only the criteria that are 
different from the signal selection criteria listed in Table ITT1 are shown. 



C. Multijet Background 

Multijet events become a background when a jet is 
misidentified as an isolated lepton or when a real lepton 
appears as a decay product of hadrons in jets, for example 
from b- or c-jets, and is sufficiently isolated. In the fol- 
lowing, such lepton-like objects are collectively referred 
to as misidentified leptons. The multijet background in 
each signal region, and in the VF+jets and tt control re- 
gions, where it is more significant, is estimated from the 
data following a matrix method similar to that employed 
in Ref. |16| . 

The multijet background from all sources (but sepa- 
rated by lepton flavor) is determined collectively. In the 
single-lepton channels, the multijet process is enhanced 
in control samples with all the signal region criteria ap- 
plied but where the lepton isolation criteria are not im- 
posed and the shower shape requirements on electrons 
are relaxed. Defining N pajSS and iVf a ji as the number of 
events in such a loose sample passing or failing the final 
lepton selection criteria, and defining -/Vr 0a i and -/V m i s id. 
as the number of real and the number of misidentified 
leptons, the following equations hold: 

cal ^misid. -^misid. ; 
Mail = (1 — £rcal)Mcal + (1 ~ Cmisid.) Mnisid. , 

where e rca i is the relative identification efficiency for real 
leptons, and e m ; s id. is the misidentification efficiency for 
misidentified leptons. Solving the equations leads to: 

, r _ Mail ~ (1/ercal - l)JV p a.ss 

JV misid. — £ misid. J Vmisid. = T~r _ -, / • 

^/fcmisid. -L/fcrcal 

The efficiency e roa i is measured from data samples of 
Z -> U decays. 

The lepton misidentification efficiency is obtained as 
follows. For electrons (muons) with px > 25 (20) GeV 
£misid. is estimated with events containing at least one 



electron (muon) satisfying the relaxed criteria, and at 
least one signal jet with px > 30 (60) GeV. In addition, 
for the electron case, i? T llss < 30 GeV is required. For 
the muon case, the event is required to contain exactly 
one muon with \do\/<Jd a > 5 where do and ad are the 
transverse impact parameter and its uncertainty, respec- 
tively, measured with respect to the primary vertex. For 
the soft-lepton channel, the sample for deriving e m i s id. 
consists of events containing a same-sign and same-flavor 
lepton pair where the leptons satisfy the relaxed isola- 
tion criteria. The selection of a lepton pair allows the 
low-px region to be studied with a large data sample. 
The same-sign requirement reduces the dominance of b- 
hadrons in the sample, providing a better mix of the dif- 
ferent mechanisms by which leptons can be misidentified. 
One of the leptons is required to fail the signal lepton 
criteria to further enhance the background; the misiden- 
tification efficiency is measured with the other lepton. 
An additional veto around the Z boson mass is applied. 
In all channels, the electron misidentification efficiency 
is evaluated separately for samples enhanced (depleted) 
in heavy-flavor contributions by requiring (vetoing) a b- 
tagged jet in the event. 

For the multi-lepton channels, the misidentification 
probabilities as determined above are applied to the num- 
ber of events where two leptons pass the loose selection 
criteria. The contribution from processes where one lep- 
ton is real and the other misidentified has been studied 
in both simulation and data, using a generalization of 
the above matrix method to two leptons. Both methods 
give similar results; the final estimate is taken from the 
simulation. 



D. Other Backgrounds 

The backgrounds from single-top, diboson and tt+ vec- 
tor boson production are estimated almost purely from 
simulation, as is the Z+jets background for the single- 
lepton channels. The background from cosmic-ray muons 
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FIG. 2. Top: m^f distribution in the W+jets (left) and ti (right) control regions for data and simulation for the single hard- 
lepton channels. Bottom: Distribution of the number of jets in the W+jets (left) and ti (right) control regions. In all plots, the 
last bin includes all overflows. The electron and muon channels are combined for ease of presentation. The "Data/SM" plots 
show the ratio between data and the total Standard Model expectation. The expectation for multijets is derived from the data. 
The remaining Standard Model expectation is entirely derived from simulation, normalized to the theoretical cross sections. 
The uncertainty band around the Standard Model expectation combines the statistical uncertainty on the simulated event 
samples with the systematic uncertainties on the jet energy scale, 6-tagging, data-driven multijet background, and luminosity. 
The systematic uncertainties are largely correlated from bin to bin. An example of the distribution for a simulated signal is 
also shown (not stacked); the signal point is chosen to be near the exclusion limit of the analysis in Ref. 



overlapping a hard-scattering event is estimated from a 
control sample with large zq, denned as the distance in 
the z direction with respect to the primary vertex, eval- 
uated at the point of closest approach of the muon to the 
primary vertex in the transverse plane. The extrapolated 
contribution to the signal region, \zq\ < 5 mm, is found 
to be negligible. 



VIII. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties have an impact on the ex- 
pected background and signal event yields in the control- 
and signal regions. These uncertainties are treated as 
nuisance parameters in a profile likelihood fit described 
in Sec. IIXI The following systematic uncertainties on 
the reconstructed objects are taken into account. The 
jet energy scale (JES) uncertainty has been determined 
from a combination of test beam, simulation and in-situ 
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FIG. 3. Top: E™ lss /m e fi distribution in the VK+jets (left) and ti (right) control regions for data and simulation for the soft- 
lepton channel. Bottom: Jet multiplicity distribution in the VK+jets (left) and ti (right) control regions. In all distributions, 
electron and muon channels are combined. The "Data/SM" plots show the ratio between data and the total Standard Model 
expectation. The expectation for multijets is derived from the data. The remaining Standard Model expectation is entirely 
derived from simulation, normalized to the theoretical cross sections. The uncertainty band around the Standard Model 
expectation combines the statistical uncertainty on the simulated event samples with the systematic uncertainties on the jet 
energy scale, 6-tagging, data-driven multijet background, and luminosity. The systematic uncertainties are largely correlated 
from bin to bin. An example of the distribution for a simulated signal is also shown (not stacked); the signal point is near the 
exclusion limit of this analysis. 



measurements from 2010 pp collision data (68[. Addi- 
tional contributions from the higher luminosity and pile- 
up in 2011 are taken into account. Uncertainties on 
the lepton identification, momentum/ energy scale and 
resolution are estimated from samples of Z — > £ + £~, 
J/tp -> £+£- and W ± l±v decays in data [HH. 
The uncertainties on the jet and lepton energies are prop- 
agated to the E™ 1SS ; an additional E™ ms uncertainty 
arising from energy deposits not associated with recon- 
structed objects is also included [74j]. Uncertainties on 



the 6-tagging efficiency are derived from dedicated data 
samples [ZaIZfl, e -g- containing muons associated with 
jets. Uncertainties on the light-flavor mis-tag rate are de- 
rived by examining tracks with negative impact param- 
eter [77| while charm mis-tag uncertainties are obtained 
from data samples tagged by reconstructing D* mesons 

Uncertainties in the matrix method for the determi- 
nation of the multijet background include the statistical 
uncertainty in the number of events available in the var- 



11 





FIG. 4. Top: mj,g c distribution in the Z+jets (left) and ti (right) control regions for data and simulation for the multi- 
lepton channels. Bottom: Distribution of the number of jets in the Z+jets (left) and ti (right) control regions; the last bin 
includes all overflows. The ee and /x/i channels are combined for Z+jets and ee, fifi and e/i channels are combined for the 
ti distributions for ease of presentation. The "Data/SM" plots show the ratio between data and the total Standard Model 
expectation. The expectation for multijets is derived from the data. The remaining Standard Model expectation is entirely 
derived from simulation, normalized to the theoretical cross sections. The uncertainty band around the Standard Model 
expectation combines the statistical uncertainty on the simulated event samples with the systematic uncertainties on the jet 
energy scale, 6-tagging, data-driven multijet background, and luminosity. The systematic uncertainties are largely correlated 
from bin to bin. An example of the distribution for a simulated signal is also shown (not stacked); the signal point is chosen 
to be near the exclusion limit of the analysis in Ref. [7^ . |73| . 



ious control samples, the difference in misidentification 
efficiency for electrons from heavy- versus light-flavored 
jets, the dependence of the misidentification efficiency 
on the jet multiplicity, and the uncertainty in the sub- 
traction of other backgrounds from the samples used to 
estimate the misidentification efficiency. 

Uncertainties from the identification efficiency for jets 
associated with the primary vertex and from the over- 
lay of pile-up in simulated events are both found to be 



negligible. 

Theoretical modeling uncertainties in the simulation 
include the following contributions. In previous versions 
of the analysis, renormalization and factorization scale 
uncertainties were estimated by varying the correspond- 
ing parameters in the ALPGEN generator by a factor 
of two, up and down from their nominal settings. Since 
these variations affect mostly the overall normalization of 
the cross sections for the samples with different values of 
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FIG. 5. Top: Distribution of the pt of the Z boson in a region enhanced in Z+jets events (ee and fifi final states combined) 
before (left) the application of any reweighting factors, and after (right) the final fit to all background control regions (described 
in Sec. IIX|I . Bottom: Distribution of _B™ 1SS in the hard-lepton VK+jets control region (electron and muon channels combined, 
lower requirement on _E™ 1SS set to 50 GeV and upper requirement removed) before (left) application of reweighting factors and 
after (right) the final fit to all background control regions. Events in the overflow bin have not been plotted. 



^Vparton, they are replaced here by a normalization of the 
individual light-parton bins to the data (see Sec. IIX|) . Ad- 
ditional generator uncertainties arise from the parameter 
that describes the jet pt threshold used in the match- 
ing (pT,min)- This uncertainty is assessed by changing 
its default value from 15 GeV to 30 GeV; the difference 
is assigned as both a positive and negative uncertainty. 
Uncertainties arising from initial- and final-state radia- 
tion are taken into account by the variation of the MLM 
matching parameter in multi-leg generators as well as by 
studying dedicated PYTHIA tunes with increased or de- 
creased radiation [79| . Fragmentation/hadronization un- 
certainties are estimated by comparing HERWIG with 
PYTHIA. In order to vary the heavy-flavor fraction, the 
cross sections for Wbb+iets and Wcc+jets in Table U 



are scaled by 1.63 ± 0.76, while VFc+jets is scaled by 
1.11 ± 0.35, based on correction factors derived from 
data [54( | . The uncertainty on Zbb+jets is taken to be 
±100%. The uncertainties on the cross sections for ti+W 
and ti+Z are taken from the NLO calculations in Refs. 




The uncertainty in the signal cross section is taken 
from an envelope of cross section predictions using dif- 
ferent PDF sets (including the as uncertainty) and fac- 
torization and renormalization scales, as described in 
Ref. (8JJ. For the simplified models, uncertainties in 
the modeling of initial-state radiation play a significant 
role for low gluino masses and for small mass differences 
in the decay cascade. These uncertainties are estimated 
by varying generator tunes in the simulation as well as 
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by generator-level studies of gg and production with an 
additional ISR jet generated in the matrix element with 
MADGRAPH5. 

The impact of these systematic uncertainties on the 
background yields and signal estimates are evaluated via 
an overall fit, described in Sec. IIXI and IXl 

IX. BACKGROUND FIT 

The background in the signal region is estimated with 
a fit based on the profile likelihood method J81]. The 
inputs to the fit are as follows: 

1. The observed numbers of events in the W+jets (or 
Z+jets in the multi-lepton channels) and tt con- 
trol regions, and the numbers expected from simu- 
lation. These are separated into 7 jet-multiplicity 
bins, ranging from 3 to 9 jets for the hard-lepton 
channel, 8 jet multiplicity bins ranging from 2 to 9 
jets for the multi-lepton channels, and 6 bins rang- 
ing from 2 to 7 jets for the soft-lepton channel. This 
information is shown in the bottom half of Fig. [2] 

tola 

2. Transfer factors (TF), derived from simulation, 
are multiplicative factors that propagate the event 
counts for VF+jets, Z+jets and tt backgrounds 
from one control region to another, or from one 
control region to the signal region. Typical values 
of the TFs from the control to the signal region are 
10~ 2 and 10~ 4 for the soft- and hard-lepton chan- 
nels, respectively. 

3. The number of multijet background events in all 
control and signal regions, as derived from the data. 

4. Expectations from simulation for the number of 
events from the minor backgrounds (single-top, di- 
boson) in all control and signal regions. 

For each analysis channel (hard-lepton, soft-lepton, 
multi-lepton) the event count in each bin of the control 
region is treated with a Poisson probability density func- 
tion. The statistical and systematic uncertainties on the 
expected yields are included in the probability density 
function as nuisance parameters, constrained to be Gaus- 
sian with a width given by the size of the uncertainty. Ap- 
proximately 150 nuisance parameters are included in the 
fit. Correlations in the nuisance parameters from bin to 
bin are taken into account where necessary. The Poisson 
probability density functions also include free parame- 
ters, for example to scale the expected contributions from 
the major backgrounds; these are described in more de- 
tail below. A likelihood is formed as the product of these 
probability density functions and the constraints on the 
nuisance parameters. Each lepton flavor (in the multi- 
lepton channel, each combination of flavors of the two 
leading leptons) is treated separately in the likelihood 
function. The free parameters and nuisance parameters 



are adjusted to maximize the likelihood. An important 
difference with respect to the analysis in Ref. [lj| is the 
increase in the number of measurements, allowing the fit 
to be constrained. This has been used in this analysis 
to constrain the nuisance parameters for the jet energy 
scale and the uncertainty in the ALPGEN scale param- 
eters from the shape information provided in the control 
regions. 

The free parameters considered in the fit are as follows: 

1. tt background: Each tt sample, broken down by 
-Wparton bin (from to 3, with the last being in- 
clusive), is scaled by a free parameter. For each 
-^Vparton bin, a common parameter is used for semi- 
leptonic and dileptonic tt samples. 

2. W/Z background: Each VF+jets and Z+jets sam- 
ple, again broken down by -/Vp ar t n bin from 2 to 
5, is scaled by a free parameter. The N pajIt on = 6 
bin for W + light-flavored jets shares its fit pa- 
rameter with -/Vp ar ton = 5. The vector boson plus 
heavy-flavor samples share the same relative nor- 
malization parameters as the light-flavor samples. 
Only iVp ar ton bins between two and five are allowed 
to float, as the lower multiplicity bins suffer from 
small numbers of events due to the jet and effective 
mass requirements. 

The backgrounds from multijets and the sub-dominant 
backgrounds from single-top and diboson production are 
allowed to float in the fit within their respective uncer- 
tainties. 

Notable nuisance parameters in the fit are: 

1. The uncertainty in the ALPGEN MLM- matching 
parameter pT,min manifests itself in the relative 
normalization of the ALPGEN -/Vp ar t on samples and 
in the jet spectra within each sample. The 
change in the event counts in the array of all con- 
trol regions, resulting from this shift in the relative 
normalization, is mapped to one parameter for both 
VF+jets and Z+jets and a separate parameter for 
tt. 

2. The uncertainty in the normalization of the 
^Vparton = 0, 1 bins for VF+jets and Z+jets, due to 
uncertainties in renormalization and factorization 
scales, is treated by one nuisance parameter. 

3. The overall normalization of the vector boson plus 
heavy flavor samples is assigned a nuisance param- 
eter reflecting the uncertainty in the cross section. 

4. The uncertainty from the fit of the p^ distribution 
is treated by assigning one nuisance parameter for 
each bin in true pr- Four equal- width bins are used 
from to 200 GeV, and one bin for p T > 200 GeV. 

5. The uncertainty due to the jet energy scale is con- 
sidered in three jet p T bins (25-40 GeV, 40-100 
GeV and > 100 GeV). The resulting uncertainty in 
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the event counts in the array of all control regions is 
mapped to one nuisance parameter for each of the 
three jet pr bins. The usage of three jet pr bins 
prevents the fit from artificially over-constraining 
the jet energy scale. 



A. Background Fit Validation 

The background fit is cross-checked in a number of val- 
idation regions, situated between the control and signal 
regions, where the results of the background fit can be 
compared to observation. These validation regions are 
not used to constrain the fit. For the single hard-lepton 
channels, one common set of validation regions, which 
receives contributions from both 3- and 4-jet channels, is 
defined as follows: 

1. The VF+jets validation region is identical to the 
TU+jets control region for the 3-jet channel ex- 
cept that the E^ lss requirement is changed to 
150 GeV < £^ iss < 250 GeV (from [40, 150] GeV). 

2. Similarly, the tt validation region is identical to the 
tt control region for the 3-jet channel except for 
the change in the E™ lss requirement to 150 GeV < 
E™ iss < 250 GeV (from [40, 150] GeV). 

3. The high transverse mass validation region is de- 
fined by m T > 100 GeV and 40 GeV < E™ iss < 
250 GeV. This region tests the validity of the back- 
ground yields from dileptonic tt events. 

For the soft-lepton channel, the validation region is based 
on the sum of the lU+jets and tt control regions but 
with the transverse mass selection changed to 80 GeV < 
m T < 100 GeV (from [40, 80] GeV). 

For the multi-lepton channels, two Z+jets validation 
regions and two tt validation regions are defined: 

1. The 2-jet Z+jets validation region is similar to the 
Z+jets control region with > 2 jets, but the leading 
two jets are required to have pr >120 GeV (instead 
of 80 GeV and 50 GeV); the fourth leading jet (if 
present) is required to have pt <50 GeV. 

2. The 4-jet Z+jets validation region is similar to the 
control region with > 4 jets but the leading jet 
Pt requirement is tightened to pt > 80 GeV (in- 
stead of 50 GeV) . 

3. The 2-jet tt validation region is similar to the 
tt control region with > 2 jets but the leading 
two jets are required to have pt > 120 GeV (in- 
stead of 80 GeV and 50 GeV); the fourth leading 
jet (if present) is required to have pt < 50 GeV. 
The E™ 1SS requirement is changed to 100 GeV < 
E^ iss < 300 GeV. 

4. The 4-jet validation region is similar to the tt con- 
trol region with > 4 jets but the leading jet pr re- 
quirement is tightened to pt > 80 GeV (instead 



of 50 GeV). The E™ lss requirement is tightened to 
80 GeV < Ef iss < 100 GeV. 

In both Z+jets validation regions the E^ lss requirement 
is tightened to 50 GeV < E™ iss < 100 GeV, and the num- 
ber of 6-tagged jets is required to be zero in order to sup- 
press the ti contamination. The selection requirements 
for the validation regions are summarized in Tables IIVI 
and [V] for the single-lepton and multi-lepton channels, 
respectively. 

The results of the fit to the control regions, as well 
as the comparison of observed versus predicted event 
counts in the validation regions, are summarized in Fig. 
|U The difference between the observed and predicted 
event counts is normalized by the total (statistical and 
systematic) uncertainty on the prediction. The agree- 
ment between predicted and observed yields is good. 



X. RESULTS AND INTERPRETATION 

The predicted background in the signal regions and the 
observed numbers of events are shown in Tables I VII and 
IVIII The data are consistent with SM expectations in all 
signal regions. 

The dominant background uncertainty comes from the 
limited number of events in the background simulation 
samples in the signal region. Uncertainties on the jet 
energy scale and the scale uncertainties for the tt back- 
ground at high jet multiplicity are also significant. In 
the soft-lepton channel, an important contribution comes 
from the evaluation of the multijet background. 

For the signal prediction, the dominant uncertainties 
at the highest excluded SUSY masses arise from the 
PDFs (30-40%) and the JES (10-20%); the former re- 
flect the uncertainty in the gluon distribution at high 
values of a;. In the simplified models with small mass 
differences typical uncertainties from ISR variations are 
approximately 30%. 

Model-independent limits on the visible cross section 
(i.e. the cross section evaluated inside a given signal re- 
gion) are derived by including the number of events ob- 
served in that region as an input to the fit and deriving 
an additional parameter, representing the non-SM signal 
strength (constrained to be non- negative) , as the output 
of the fit. Potential signal contamination in the con- 
trol regions is ignored. Limits on the number of non-SM 
events in the signal region, derived using the CL S 82 1 
prescription, are divided by the integrated luminosity to 
obtain the constraints on the visible cross section. The 
limits at 95% confidence level (CL) are shown in Table 

Em] 

For excluding specific models of new physics, the fit in 
the signal regions proceeds in the same way except that 
in this case the signal contamination in control regions 
is treated by providing transfer factors from the signal 
regions to the control regions as further input to the fit. 
In addition, the likelihood fit makes use of the rn"^ shape 
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hard-lepton 




soft-lepton 




W VR 


ti VR 


High-m T VR 




N ie t 


> 3 


> 3 


> 3 


> 2 


(GeV) 


> 80, 25, 25 


> 80, 25, 25 


> 80, 25, 25 


> 130,25 


iVjet (6-tagged) 





> 1 






E™ iss (GeV) 


[150,250] 


[150,250] 


[40,250] 


[180,250] 


m T (GeV) 


[40,80] 


[40,80] 


> 100 


[80,100] 


m£f (GeV) 


> 500 


> 500 


> 500 





TABLE IV. Overview of the selection criteria for the background validation regions (VR) for the single-lepton channels. Only 
the criteria that are different from the signal selection criteria listed in Table HI] are shown. 



multi-lepton 2-jet multi-lepton 4-jet 

Z VR ti VR Z VR ti VR 



ATjet 


> 2 


> 2 


> 4 


> 4 


(GeV) 


> 120, 120 


> 120, 120 


> 80,50,50,50 


> 80,50,50,50 


Njet (fe-tagged) 




> 1 




> 1 


E!f iss (GeV) 


[50,100] 


[100,300] 


[50,100] 


[80,100] 


mu (GeV) 


[81,101] 


< 81 or > 101 


[81,101] 


< 81 or > 101 



TABLE V. Overview of the selection criteria for the background validation regions (VR) for the multi-lepton channels. Only 
the criteria that are different from the signal selection criteria listed in Table ITT1 are shown. For the 2-jet validation regions, the 
fourth leading jet (if present) is required to have pt < 50 GeV. 



information (E™ lss / m e s for the soft-lepton channel) in 
the signal region as a further discriminant. Examples 
of these distributions are shown in Fig. [7] (the figure 
shows the distributions summed over lepton flavors, while 
the fit treats each lepton flavor channel independently). 
The likelihood is extended to include bin-by-bin m"^ c or 
E™ lss /m c H information by dividing the signal region into 
several bins of m 1 ^ or E^ lss /m c s. 

The ten statistically independent hard-lepton and 
multi-lepton channels are combined to set limits in the 
MSUGRA/CMSSM model. For the minimal GMSB 
model, only the multi-lepton channels are used. The soft- 
lepton channels are used together with the hard-lepton 
and multi-lepton channels to set limits in the one- and 
two-step simplified models. 

The limit in the plane of versus mo in the 

MSUGRA/CMSSM model is shown in Fig. E The 
band around the expected limit includes all uncertainties 
except theoretical uncertainties on the signal prediction 
while the band on the observed limit indicates the sen- 
sitivity to the theoretical uncertainties on the signal. A 
large improvement in exclusion coverage over the previ- 
ous analysis [iH ] can be seen. The simultaneous fit to 
the ten signal regions and the inclusion of the shapes 
of the Wg'jj distributions increase the expected reach in 
m-iM and mo by about 100 GeV, approximately uniformly 
across the plane. Along the line of equal masses between 
squarks and gluinos in the MSUGRA/CMSSM model, 



masses below approximately 1200 GeV are excluded at 
95% CL. 

For the minimal GMSB model, the limit in the plane of 
tan /3 versus A is shown in Fig. [SJ The exclusion reach is 
dominated by the dilepton plus two jets channel. Values 
of A below about 50 TeV are excluded at 95% CL for 
tan/3 < 45, improving on previous constraints. 

Exclusion limits in the one-step simplified models are 
shown in Fig. [TUJ The figures also show the cross section 
excluded at 95% CL. The exclusion limits in the two- 
step simplified models are shown in Fig. [11] for gluino 
pair production and Fig. [T2]for squark pair production. 
Simplified models with varying chargino mass and two- 
step simplified models are considered here for the first 
time in leptonic SUSY searches. For both one- and two- 
step models, for the case of low LSP masses, gluinos with 
masses below approximately 900-1000 GeV and squarks 
with masses below approximately 500-600 GeV are ex- 
cluded. Squark limits are considerably weaker, primarily 
due to the lower production cross section. Furthermore in 
the one-step model, gluinos with mass below 550 GeV are 
excluded for essentially all values of the LSP mass if the 
latter is more than 30 GeV smaller than the gluino mass. 
Care has to be taken when interpreting the simplified 
model limit in the context of a MSSM scenario, where 
the mass of the sneutrino is lighter than the mass of the 
left-handed slepton, as this can lead to modification of 
the lepton momenta. 
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FIG. 6. Summary of the fit results in the control regions (left) and validation regions (right). The difference between the 
observed and predicted number of events, divided by the total (statistical and systematic) uncertainty on the prediction, is 
shown for each control and validation region. 
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Single-lepton 




Electron 








Muon 




Number of events 


3- jet 


4-jet 


soft lepton 


3- jet 


4-jet 


soft lepton 


Observed 


2 


4 




11 


1 


2 


14 


Fitted bkg 


2.3 ± 0.9 


3.5 ± 0.9 


14 


.0 ± 3.3 


2.6 ± 0.8 


1.5 ± 0.3 


19 ± 5 


Fitted top 


0.4 ±0.2 


2.3 ±0.6 


3. 


8 ±0.6 


0.5 ±0.2 


1.3 ±0.3 


3.8 ±0.8 


Fitted W/Z+jets 


1.5 ±0.6 


0.9 ±0.2 


5. 


8± 1.0 


2.0 ±0.6 


0.2 ±0.1 


11.4 ±2.3 


Fitted other bkg 


0.0 ±0.0 


0+ ' 3 


0. 


6 ±0.1 


0.1 ±0.1 


0.0 ±0.0 


0.2 ±0.1 


Fitted multijet 


0.3 ±0.4 


0.3 ±0.4 


3. 


8 ±2.5 


0.0 ±0.0 


0.0 ±0.0 


3.6 ±2.5 


MC exp. SM 


2.7 


5.3 




14.2 


2.8 


2.4 


18.0 


MC exp. top 


0.9 


3.1 




4.3 


0.6 


2.0 


3.8 


MC exp. W/Z+jets 


1.5 


1.3 




5.5 


2.0 


0.3 


10.5 


MC exp. other bkg 


0.0 


0.5 




0.5 


0.2 


0.1 


0.1 


Data-driven multijet 


0.3 


0.3 




3.8 


0.0 


0.0 


3.6 



TABLE VI. The observed numbers of events in the single-lepton signal regions, and the background expectations from 
the fit. The inputs to the fit are also shown; these consist of the data-driven multijet background estimate and the nominal 
expectations from simulation (MC), normalized to theoretical cross sections. The errors shown are the statistical plus systematic 
uncertainties. 



Multi-lepton 




2-jets 






4-jets 




Number of events 


ee 


HH 


e/i 


ee 




e/i 


Observed 








1 


8 


12 


18 


Fitted bkg 


0.3 ±0.2 


0.4 ± 0.2 


0.7 ±0.2 


9.1 ± 1.5 


11.7 ± 1.7 


21 ± 3 


Fitted top 


0.1 ±0.1 


0.2 ±0.1 


0.6 ±0.2 


9.1 ± 1.4 


11.1 ± 1.7 


20 ±3 


Fitted W/Z+jets 


0.1 ±0.1 


0.1 ±0.0 


0.0 ±0.0 


0.0 ±0.0 


0.2 ±0.1 


0.4 ±0.1 


Fitted other bkg 


0.1 ±0.1 


0.1 ±0.0 


0.1 ±0.0 


0.0 ±0.0 


0.4 ± 0.1 


0.6 ±0.1 


Fitted multijet 


0.0 ±0.0 


0.0 ±0.0 


0.0 ±0.0 


0.0 ±0.2 


0.0 ±0.0 


0.0 ±0.0 


MC exp. SM 


0.3 


0.5 


0.9 


11.4 


14.7 


27.1 


MC exp. top 


0.2 


0.3 


0.7 


11.1 


13.9 


26.0 


MC exp. W/Z+jets 


0.1 


0.1 


0.1 


0.1 


0.3 


0.4 


MC exp. other bkg 


0.1 


0.1 


0.1 


0.2 


0.5 


0.7 


Data-driven multijet 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 



TABLE VII. The observed numbers of events in the multi-lepton signal regions, and the background expectations from 
the fit. The inputs to the fit are also shown; these consist of the data-driven multijet background estimate and the nominal 
expectations from simulation (MC), normalized to theoretical cross sections. The errors shown are the statistical plus systematic 
uncertainties. 
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Signal channel 




C.95 
"Jobs 


o95 




hard electron, 3-jet 


0.94 


4.4 


4 O+2.0 


0.54 


hard muon, 3-jet 


0.75 


3.6 


4.2t 2 ;? 


0.27 


hard electron, 4-jet 


1.22 


5.8 




0.63 


hard muon, 4-jet 


0.95 


4.5 


o Q+1.3 
d.8_ .7 


0.75 


soft electron 


1.82 


8.6 




0.28 


soft muon 


1.92 


9.0 


12.5+J1 


0.21 


multi-lepton, ee, 2-jet 


0.71 


3.3 


3.5 ±0.1 


0.48 


multi-lepton, ////, 2-jet 


0.76 


3.6 


3.5 ±0.1 


0.46 


multi-lepton, e/j,, 2-jet 


0.83 


3.9 


3 6+ ' 6 


0.85 


multi-lepton, ee, 4-jet 


1.53 


7.2 


fy 7+3.2 

'■'-2.1 


0.39 


multi-lepton, fifx, 4-jet 


1.93 


9.1 




0.55 


multi-lepton, e/i, 4-jet 


2.14 


10.1 


n.5±|;l 


0.28 



TABLE VIII. Left to right: 95% CL upper limits on the visible cross section ((ecr)^) m the various signal regions, and on 
the number of signal events {S'ihs )• The third column (S'cxp) shows the 95% CL upper limit on the number of signal events, 
given the expected number (and ±la uncertainty on the expectation) of background events. The last column indicates the 
CLb value, i.e. the observed confidence level for the background-only hypothesis. 
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FIG. 7. Top and middle: Distribution of mjjj in the signal regions after all selection requirements except for that on the 
inclusive effective mass. Top left: hard-lepton, 3-jet selection. Top right: hard-lepton, 4-jet selection. Middle left: multi- 
lepton, 2-jet selection. Middle right: multi-lepton, 4-jet selection. The last m 1 ^ bin includes all overflows. The lowest mjjfj? 
bins are affected by the minimum pt requirements on jets and _E™ 1SS . Bottom: The E™ laa /m e ff distribution in the soft-lepton 
signal region after all selection requirements except for that on E™ 1SB /m g. In all plots the different lepton flavors have 
been combined for ease of presentation. The "Data/SM" plots show the ratio between data and the total Standard Model 
expectation. The Standard Model expectation shown here is the input to the final fit, and is derived from simulation only, 
normalized to the theoretical cross sections. The uncertainty band around the Standard Model expectation combines the 
statistical uncertainty on the simulated event samples with the systematic uncertainties on the jet energy scale, 6-tagging, 
data-driven multijet background, and luminosity. The systematic uncertainties are largely correlated from bin to bin. An 
example of the distribution for a simulated signal is also shown (not stacked) ; the signal point is chosen to be near the exclusion 
limit of the analysis in Ref. [la ]. 
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mSUGRA/CMSSM: tanp = 1 0, A = 0, n > 
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FIG. 8. Expected and observed 95% CL exclusion limits in the MSUGRA/CMSSM model with tan/3 = 10, A Q = and 
the sign of /i taken to be positive. The results are obtained by combining ten signal regions from the hard single-lepton and 
multi-lepton channels. The band around the median expected limit shows the ±lcr variations, including all uncertainties except 
theoretical uncertainties on the signal. The dotted lines around the observed limit indicate the sensitivity to ±l(j variations on 
these theoretical uncertainties. The dashed grid shows contours of constant squark (curved lines) and gluino (nearly horizontal 
lines) masses. The previous limit from ATLAS fl6l ] and the results from the LEP experiments [83l ] are also shown. 




FIG. 9. Expected and observed 95% CL exclusion limits in the minimal GMSB model, combining six signal regions from the 
multi-lepton channels. The band around the median expected limit shows the ±1<t variations, including all uncertainties except 
theoretical uncertainties on the signal. The dotted lines around the observed limit indicate the sensitivity to ±lcr variations on 
these theoretical uncertainties. The different next-to-lightest-SUSY particle (NLSP) regions are indicated. The coNLSP region 
denotes the region where f\ and £r are nearly mass degenerate. Previous OPAL and ATLAS limits in this model can be found 
in Refs. [84[ and [7^. |73|. respectively. Limits derived from the LEP slepton mass limits are also shown. 
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FIG. 10. Excluded regions at 95% confidence level in the parameter space of one-step simplified models. Top row: gluino pair 
production with g — > qq'xf — > qq'W ± x1. Bottom row: squark pair production with qi, — > q'xi ~> Q , 'W /± Xi- ^ n tne left column, 
the chargino mass is set to be halfway between gluino (top) or squark (bottom) and LSP masses. In the right column, the LSP 
mass is fixed at 60 GeV and the masses of the chargino and gluino (top) or squark (bottom) are varied. The band around the 
median expected limit shows the ±ltr variations, including all uncertainties except theoretical uncertainties on the signal. The 
dotted lines around the observed limit indicate the sensitivity to ±la variations on these theoretical uncertainties. The plots 
are from the combination of the hard and soft single-lepton channels. The numbers indicate the excluded cross section in fb. 
A smaller excluded cross section implies a more stringent limit. 



22 



g-g decays via sleptonsigg-xjqqqx^-Hiqqqllvvjc^ 



I i i i i I i i i 

ATLAS 

J 



1 1 1 1 1 1 1 



1 1 1 1 1 1 



Ldt = 4.7fb ,Vs=7TeV 



__ Observed limit (+1 oJJJJ* 
Expected limit (±1 o exp ) 

All limits at 95% CL 



^- 1200 

in > — i 

— S E 1000 

- w 



800 



g-g decays via WZ: gg-Hjqqqx^-^qqqWZWZx^ 




I i i i i I i i i i 

ATLAS 



300 400 500 600 700 800 900 1000 1100 1200 

m § [GeV] 



600 



400 



200 



I " " I 



I 1 1 1 1 I 



Ldt = 4.7fb ,Vs=7TeV 



■19. 



1565 



, Observed limit (+1 <r 



1279 675 

■ 13629 1 476 460 

1231 719 

1331 395 — 



All limits at 95% CL , ''23526 

^2 .-■' 1102 756 



313 58 




m> Hh iff, m 1 1 1 1 1 i 6 p 1 1 1 1 1 alL 1 1 iM i'f 1 1 1 



300 400 500 600 700 800 900 1000 1100 1200 

m 3 [GeV] 



FIG. 11. Excluded regions at 95% confidence level in the parameter space of two-step simplified models with gluino pair 
production. Left: both gluinos decay via g — > qq'xi — > qq'l^UL — > qq'£ "Xi or 9 ~> Qq'Xi ~ ► w'^i ~~ > qq'vi Xi- Right: 
both gluinos decay via g — > qqXi ~ * q(fW X2 ~~ ^ Z^xl- The band around the median expected limit shows the 

±1(7 variations, including all uncertainties except theoretical uncertainties on the signal. The dotted lines around the observed 
limit indicate the sensitivity to ±lrr variations on these theoretical uncertainties. The plots are dominated by the multi-lepton 
channels. The numbers indicate the excluded cross section in fb. A smaller excluded cross section implies a more stringent 
limit. 



23 



> 

CD 

CD 



1200 



q-q decays via sleptons: qq->qqx + X + ->qq |lvv X X 



E 1000 



800 



600 



400 



200 k 



I 1 1 1 1 I 1 1 1 1 I 1 1 11 I 1 1 1 1 I 

ATLAS 

J Ldt = 4.7fb"',Vs=7TeV 



Expected limit (±1 o exp ) 

All limits at 95% CL 



0-0 

1*1 



1 1 1 1 1 1 



150 56 




> 



300 400 500 600 700 800 900 1000 1100 1200 

m~ [GeV] 



1200 



q-q decays via sleptons: qq->qqx x 



1 E 1000 



800 



600 



400 



200 



I 1 1 1 1 I 1 1 1 1 I 1 1 11 I 1 1 1 1 I 

ATLAS 

J Ldt = 4.7fb" , ,tfs=7TeV 

Observed limit (+1 of h ^) 

Expected limit (±1 o exp ) 

All limits at 95% CL 



*qq"ivx°x°,qqivvvx°x° 



i 1 1 1 1 1 




300 400 500 600 700 800 900 1000 1100 1200 

m~ [GeV] 



> 

a> 
(3 



1200 



q-q decays via WZ: qq^qqx*x*^qqWZWZx°x° 



E 1000 



800 



600 



400 - 



200 tj ; 



! | I I I I | I I I I | 

ATLAS 

J Ldt»4.7fb' 1 ,Vs=7TeV 



I " " I 



I " " I 



. Observed limit (+1 a™) 
■ Expected limit (±1 o exc ) 



All limits al 95% CL 




300 400 500 600 700 800 900 1000 1100 1200 

m~ [GeV] 



qvv L 



qwXi 



q XT ~^ q'l v — > q'l^fXi an d the other 



FIG. 12. Excluded regions at 95% confidence level in the parameter space of two-step simplified models with squark pair 
production. Top left: both squarks decay via qL — > q'xt ~* q'l vl — > q'l vXi or Ql — > q'xt ~ * q'l v 
Top right: one squark decays via qL — > q'xt ~* q'l Vl — > q'l VX^ or qL ' J: 

squark decays via -> q%% — > q£ l^ — > ql^l^Xi or 

q L q'xt W^xl : ' ~ " 
uncertainties except theoretical uncertainties on the signal. The dotted lines around the observed limit indicate the sensitivity 
to ±1ct variations on these theoretical uncertainties. The plots are dominated by the multi-lepton channels. The numbers 
indicate the excluded cross section in fb. A smaller excluded cross section implies a more stringent limit. 



Bottom row: both squarks decay via 



' - >• W^^x® — > W^^Z^Xi- The band around the median expected limit shows the ±1<j variations, including all 



24 



XI. CONCLUSION 

A new search with the ATLAS detector for SUSY in 
final states containing jets, one or more isolated lep- 
tons (electron or muon) and i?™ ISS has been presented. 
Data from the full 2011 data-taking period, correspond- 
ing to an integrated luminosity of 4.7 fb _1 , have been 
analyzed. Single- and multi-lepton channels are treated 
in one analysis. A signal region with a soft lepton and 
soft jets has been introduced to increase the sensitivity 
to SUSY decay spectra involving small mass differences 
("compressed SUSY"), where the sensitivity is improved 
by a factor of 10-30 compared to the hard-lepton chan- 
nel. A simultaneous fit is performed to the event yield in 
multiple signal and control regions and to the shapes of 
distributions in those regions. 

Observations are in good agreement with SM expecta- 
tions and constraints have been set on the visible cross 
section for new physics processes. Exclusion limits have 
also been extended for the MSUGRA/CMSSM and min- 
imal GMSB models as well as for a number of simpli- 
fied models. In MSUGRA/CMSSM, squark and gluino 
masses below approximately 1200 GeV are excluded at 
95% CL (for equal squark and gluino masses). In mini- 
mal GMSB, values of A below about 50 TeV are excluded 
for tan/3 < 45. 

In one-step simplified models (with the chargino mass 
halfway between the masses of the gluino/squark and 
LSP) gluinos are excluded for masses below approxi- 
mately 900 GeV for low values of the LSP mass. Gluinos 
with mass below 550 GeV are excluded for essentially 
all values of the LSP mass if the latter is more than 
30 GeV smaller than the mass of the gluino. In the one- 
step simplified model with a fixed LSP mass and varying 
chargino and gluino (squark) masses, gluinos below ap- 
proximately 950 GeV are excluded for a wide range of 
chargino masses; squarks are excluded below 500 GeV, 
albeit for a narrower range of chargino masses. 

A variety of two-step simplified models have been con- 
sidered. Limits on gluino masses range from about 



900 GeV to 1000 GeV, while squark mass limits range 
from about 500 GeV to 600 GeV, all for low LSP masses. 

These results improve significantly on previous con- 
straints. 
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